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Abstract 

We consider the static Maxwell system with an axially symmetric 
dielectric permittivity and construct complete systems of its solutions 
which can be used for analytic and numerical solution of corresponding 
boundary value problems. 



1 Introduction 

Consider the static Maxwell system 

div(eE) = 0, rotE = (1) 

where we suppose that e is a function of the cylindrical radial variable r = 
\/ x\ + x\: e = e(r). Two important situations are usually studied: the 
meridional field and the transverse field. 
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The first case is characterized by the condition that the vector E is inde- 
pendent of the angular coordinate 9 and the component Eg of the vector E in 
cylindrical coordinates vanishes identically The vector of such field belongs 
to a plane containing the axis X3 and depends only on the distance r to this 
axis as well as on the coordinate X3. The field then is completely described 
by a two-component vector-function in the plane (r, X3). 

The second case is characterized by the condition that the vector E is in- 
dependent of X3 and the component £"3 is identically zero. The vector of such 
field belongs to a plane perpendicular to the axis X3 and the corresponding 
model reduces to a two-component vector-function in the plane (xi,X2). 

In the present work in both cases we construct a complete system of 
solutions of the corresponding model. We use the fact that in both cases the 
system ([1]) reduces to a system describing so-called p-analytic functions pQ , 

m, h, ma, nu, m\,m 

1 1 

U X = ~Vy, Uy = ~~ V X- (2 J 

In the first case the function p is a function of one Cartesian variable x 
meanwhile in the second it is a function of r = yVx 2 + y 2 . In both cases 
we construct an infinite system of so-called formal powers [2], [5]. This is a 
complete system of exact solutions of equations (jSJ) generalizing the system 
of usual complex powers (z — z ) n , n = 0, 1, 2, . . .. Locally, near the center z 
the formal powers behave asymptotically like powers. Nevertheless in general 
their behaviour can be arbitrarily different from that of powers but with a 
guarantee of their completeness in the sense that any solution of the consid- 
ered equations can be represented as a uniformly convergent series of formal 
powers. The general theory of formal powers was developed by L. Bers [2] as 
a part of his pseudoanalytic function theory. Its application was restricted 
by the fact that only for a quite limited class of pseudoanalytic functions the 
explicit construction of formal powers was possible. L. Bers' results allow 
us to construct a complete system of formal powers in the meridional case. 
Nevertheless they are not applicable to the model arising from the trans- 
verse case. In the recent works [12] and [15] the class of solvable in this sense 
systems © was substantially extended. In the present work we use these 
results for solving the static Maxwell system in the transverse case. This 
combination of the relation between the static Maxwell system ([I]) and the 
system (JSJ together with the classical results of L. Bers on pseudoanalytic 
formal powers and new developments in [12] and [TS] allow us to obtain a 
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general solution of the static Maxwell system in the axially symmetric case 
in the sense that we construct a complete system of its solutions for both the 
meridional and the transverse fields. 



2 Reduction of the static Maxwell system to 
^-analytic functions 

2.1 The meridional case 

Introducing the cylindrical coordinates and making the assumptions that 
E is independent of the angular variable 9 and that the component Eg is 
identically zero we obtain that ([1]) can be written as follows 

^_^ = 1 d(reE r ) | d(eE 3 ) _ Q 
dx^, dr r dr dx 3 

Denote x = r, y = a; 3 , u = E 3 and v = reE r . Then the system takes the 
form 

1 1 

^x 7 \ V> } 

xs(x) xe(x) 
where the subindices denote the derivatives with respect to the corresponding 
variables. Thus, in the case of a meridional field the vector E is completely 
described by an xe(x)-analytic function uj — u + iv. 



2.2 The transverse case 

We assume that E is independent of the longitudinal variable £3 and E3 = 0. 
Then from ([ID we have that the vector [Ex, -E^) 7 is the gradient of a function 
u = u(xi 1 X2) which satisfies the two-dimensional equation 

div(Wu) = 0. (3) 

Denote x — xi, y = X2, z — x + iy and consider the system 

Ux = -v y , u y = ~-v x - (4) 

It is easy to see that if the function u = u + iv is its solution then u is 
a solution of fl3]), and vice versa [HJ, if u is a solution of ([3]) in a simply 
connected domain Q then choosing 

v = A(ieu-z), (5) 
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where 

A [$] (x, y) = 2 ^ Re $cfe + Im Qdy^J + c, (6) 

c is an arbitrary real constant, T is an arbitrary rectifiable curve in Q leading 
from (xo,yo) to (x, y) we obtain that u = u + iv is a solution of (@J. Here 
the subindex z means the application of the operator dz = \{d x + id y ). Note 
that due to the fact that u is a solution of (El) the function $ = ieu-z satisfies 
the condition d y Re $ — Im $ = and hence the integral in ([6]) is path- 
independent. For a convex domain then expression ([6]) can be written as 
follows 

A [$] 0, y) = 2 ( f Re $(77, y)^ + / Im $(x , £)d£ J + c. 
Note that w is a solution of the equation 

div(iv^) = 0. 

Thus, equation (J3]) (and hence the system ([TJ in the case under consideration) 
is equivalent to the system (jlj) in the sense that if u = u + iv is a solution of 
(j3J) then its real part w is a solution of ([2]) and vice versa, if u is a solution 
of ([3]) then u} = u + iv, where v is constructed according to (jSJ) is a solution 
of®. 

We reduced both considered cases the meridional and the transverse to 
the system describing p-analytic functions. In the first case p = xe(x) is a 
function of one Cartesian variable and in the second p = e(r), r = ^Jx 2 + y 2 . 
As we show below in both cases we are able to construct explicitly a complete 
system of formal powers and hence a complete system of exact solutions of 
the corresponding Maxwell system. Let us notice that equation ([3]) with e 
being a function of the variable r was considered in the recent work [6] whith 
applications to electrical impedance tomography. The algorithm proposed 
in that work implies numerical solution of a number of ordinary differential 
equations arising after a standard separation of variables. Our construction 
of a complete system of solutions of (j3J) is based on essentially different ideas 
and does not require solving numerically any differential equation. 
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3 ^-analytic functions and formal powers 



3.1 The main Vekua equation 

Consider the system describing p-analytic functions 

u x = ~v y , u y = — v x , (7) 
p p 

where we suppose that p is a positive and continuously differentiable function 
of x and y. Together with tis system we consider the following Vekua equation 
which due to its importance in relation to second-order elliptic equations of 
mathematical physics is called [H], [15] the main Vekua equation 

Wt = f -jW, (8) 

where / = yjp. The function uj = u + iv is a solution of ([7j) iff [15] W = 
uf + iv/fisa solution of (jSJ). 

In [T5] there was proposed a method for explicit construction of the system 
of formal powers corresponding to the main Vekua equation under a quite 
general condition on /. Here we briefly describe the method for which we 
need first to recall some basic definitions from L. Bers' theory of formal 
powers [2]. 

Let F and G be a couple of solutions of a Vekua equation 

= a {F , G) W + b {F , G) W in Q (9) 

such that lm(FG) > 0. Then (F, G) is said to be a generating pair corre- 
sponding to (jSJ). The complex functions a(p t cf) and &(f,g) are called charac- 
teristic coefficients of the pair (F, G) and it can be seen that 

FG~z — F%G FG-z — F%G 

«W = ~ fg-FG' 6(F ' G) = FG-FG- 

Together with these characteristic coefficients another pair of characteristic 
coefficients is introduced in relation to the notion of the (F, G)-derivative: 

FG Z — F Z G FG Z — F Z G 

A(F ' G) = -FG-FG' B{F > G) = FG-FG' 
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where the z means the application of the operator d z = ~ — i-J^J . As in 

the present work we do not use explicitly the notion of the (F, G) -derivative, 
we refer the interested reader to [2] for its definition and properties. However 
we do need the concept of characteristic coefficients for defining the following 
important object. 

Definition 1 Let (F,G) and (Fi,Gi) - be two generating pairs in Q corre- 
sponding to the Vekua equations with coefficients 0(f,g), °{f,g) o,nd ol{f x ,g{) 
and 6(Fi,Gi) respectively. Then (Fi,G\) is called successor of (F,G) and 
(F, G) is called predecessor of (Fi, G%) if 

a {Fx,G{] = a (F,G) and fyFi,Gi) = —B(f,g)- 

Definition 2 A sequence of generating pairs {(F m , G m )}, m = 0, ±1, ±2, . . . 
, is called a generating sequence if (F m+1 ,G m+ i) is a successor of (F m ,G m ). 
If(F ,G ) = (F,G), we say that (F,G) is embedded in {(F m ,G m )}. 

For any generating pair (F, G) the corresponding (F, G)-integral is defined 
as follows 

J wd { F,G)C = F{z) Re J G*wd( + G{z) Re J F*wd( 

where T is a rectifiable curve leading from z to z and (F*, G*) is an adjoint 
generating pair defined by the equations 

2F ^ 2G 

F* = = — G = — = — . 

FG-FG FG-FG 

If w is an (Fi, G\) - pseudoanalytic function (i.e., it is a solution of the Vekua 
equation with the coefficients a^Gi) an d &(Fi,Gi)) then its (F, G)-integral is 
path- independent . 

Now we are ready to introduce the definition of formal powers. 

Definition 3 The formal power Zm\a, z ; z) with center at z £ ^, coeffi- 
cient a and exponent is defined as the linear combination of the genera- 
tors F m , G m with real constant coefficients X, fi chosen so that XF m (zo) + 
fiG m (z ) = a. The formal powers with exponents n = 1, 2, . . . are defined by 
the recursion formula 

Z£\a,za;z) = n j (a, z ; ()d {FmiGm) (. (10) 
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This definition implies the following properties. 

1. Zm\a, z ; z) is an (F m , G m )-pseudoanalytic function of z. 

2. If a' and a" are real constants, then Zm (a'+ia", Zq; z) = a!Zm(\, Zq; z) + 
a"Zjn\i, z ;z). 

3. The asymptotic formulas 

Z%\a,z ;z)~a(z-z ) n , z - z (11) 

hold. 

Writing Z^ n '(a, zq\ z) we indicate that the formal power corresponds to 
the generating pair (F, G). 

The definition of formal powers shows us that in order to obtain Z^ n > (a, zq \ z) 
we need to have first the formal power z[ n X \a, zq; z) for which it is necessary 
to calculate Z^~ 2 \a^ z ; z) and so on. Thus, the problem of construction of 
formal powers of any order for a given generating pair (F, G) reduces to the 
construction of a corresponding generating sequence. Then definition [3] gives 
us a simple algorithm for constructing the formal powers. In other words, 
one needs a pair of exact solutions for each of the infinite number of Vekua 
equations corresponding to a generating sequence. 

In the next subsection we show how this seemingly difficult task can be 
accomplished in a quite general situation. Meanwhile here we recall some well 
known results in order to explain that the system of formal powers in fact 
represents a complete system of solutions of a corresponding Vekua equation. 
First of all, let us notice that due to the property 2 of formal powers for every 
n (and for a fixed Zq) it is sufficient to construct only two formal powers: 
■,Zq]z) and Z^ n \i, zq; z), then for any coefficient a the corresponding 
formal power Z^ n \a, zq\ z) is a linear combination of the former two. 

An expression of the form X]^=o Z^^n, zq\ z) is called a formal polyno- 
mial. Under the conditions imposed in this work on the function e and on 
the domain of interest Q (see section 4 and for more details [13] ) the fol- 
lowing Runge-type theorem is valid where following [2] we say that a series 
converges normally in a domain Q if it converges uniformly on every bounded 
closed sub domain of Q. 

Theorem 4 J&jj A pseudoanalytic function defined in a simply connected do- 
main can be expanded into a normally convergent series of formal polynomi- 
als. 
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In other words a pseudoanalytic function can be represented as an infinite 
linear combination of the functions 

{Z^(l,z ;z), ZW(i,z ;z)}™ =0 . 

Moreover, if we know that a pseudoanalytic function W satisfies the Holder 
condition on the boundary of a domain of interest Q (a common require- 
ment when a boundary value problem is considered) then, e.g., the following 
estimate in the C(fi)-norm is available. 



Theorem 5 [Tbj Let W be a pseudoanalytic function in a domain Q bounded 
by a Jordan curve and satisfy the Holder condition on dQ with the exponent 
a (0 < a < 1). Then for any e > and any natural n there exists a 
pseudopolynomial of order n satisfying the inequality 

\W{z)-P n (z)\<-^- for any z eU 

n a-e 

where the constant does not depend on n, but only on e. 

These and other results on interpolation and on the degree of approxima- 
tion by pseudopolynomials which can be found in the vast bibliography ded- 
icated to pseudoanalytic function theory (see, e.g., [7j, [9]) show us that the 
system of formal powers is as good for solving corresponding boundary value 
problems as is the system of usual complex powers (z — z ) n , n = 0, 1, 2, . . .. 
The real (or imaginary) parts of {(z — zo) n }™ =0 are harmonic polynomials 
successfully applied to the numerical solution of boundary value problems 
for the Laplace equation. In a similar way the real parts of formal powers 
|^( n )(l ; z ] z), Z( n \i, z ; z)} n _ Q corresponding to the main Vekua equation 
(JHJ) (where / = \fe) can be used for the numerical solution of boundary value 
problems for the conductivity equation ([3]) because as was shown in [T3] the 
system of functions 



\-^ReZ {n \l,z ;z), ^ReZ (n) ( 



i, zq; z / 



n=0 



is complete in the space of solutions of ([3]) in the sense of theorems H] and [5j 
A formal power (a, zq] z) related to the Vekua equation (JHD corre- 
sponds to a formal power if Z (jl \a ) z ; z) (we use the notation of L. Bers) 
related to the system (JTj) in the following way 

*Z {n) (a,z -z) = jReZ in \a,z -z)+iflmZ {n \a,z Q -z). 
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As we will see in the meridional case it is convenient to work directly with 
formal powers *Z( n \a, z ; z). Any solution to = u + iv of the system ([7]) can 
be expanded into a normally convergent series of real linear combinations of 
the complex functions {*Z( ra )(l, z ; z), *Z( n \i, Zq; z)}. 



3.2 Construction of generating sequences 

In [15] the following result was obtained. 

Theorem 6 Let F = U(u)V(v) and G = v ^ v ^ where U and V are ar- 
bitrary differentiable nonvanishing real valued functions, $ = u + iv is an 
analytic function of the variable z = x + iy in Q such that $ z is bounded 
and has no zeros in Q. Then the generating pair (F,G) is embedded in the 
generating sequence (F m , G m ), m — 0, ±1, ±2, . . .in Q defined as follows 

F m = ($ z ) m F and G m = (<& z ) m G for even m 

and 

F m = ^-F and G m = (^ z ) m U 2 G for odd m. 

This theorem opens the way for construction of generating sequences and 
consequently of formal powers in a quite general situation (see [15]) and in 
particular in both cases considered in the present work. In the meridional 
case the theorem reduces to the result of L. Bers [2| which we use in the 
next subsection while in the transverse case this and other classical results 
are insufficient for constructing formal powers explicitly and theorem [6] is 
indispensable. 



4 Construction of formal powers 

4.1 Formal powers in the meridional case 

As was shown in subsection 2.1 in the meridional case the Maxwell system 
reduces to the following couple of equations 



xe(x) xe(x) 
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which is equivalent to the system considered in [2J N18.1] 

a(x)(f) x = T(y)ip y , a(x)(f) y = -r(y)i/j x . 

Taking a(x) = xe(x) and r = 1 we can use the elegant formulas for the 
generating powers obtained by L. Bers. Let 

X^(x ,x)=X^(x ,x) = l 
and for n = 1, 2, ...denote 

X 

X (n \x ,x) = nf X (n - 1} (x ,t)^-cit for odd n 

J t£[t) 

X 

=„/xC->( Io , W for eve. 



X'(j 



X 

X {n \x ,x) = nj X^V{x ,t)te(t)dt for odd 



n 



n 



X 



I W (x , x) = n [ X {n ~ l) {x , t) -^—dt for 

J t£(t) 



even n 



Then the formal powers in the meridional case are given by the expressions 



,Z<»V + ia", z ; z) = a'J2 (f)X^ k \ k (y - y ) k 

k=0 ^ ' 

+ '«"E(^ M 'l!/-!/o) 

k=0 ^ ' 



for odd n 



and 
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,Z<»V + ia", z ; z) = a'£ Q X^ fc (?/ - y ) k 

+ ™"f;( n ^x {n - k) i k (y-yo) k for 

I— n \ / 



even n. 



4.2 Formal powers in the transverse case 

As was shown in subsection 2.2 the Maxwell system ([1]) in the transverse case 
reduces to the system 



1 1 

II r Zly V X 



where e is a positive differentiable function of r = a/x 2 + y 2 . This system 
describing £-analytic functions is equivalent to the main Vekua equation (151) 
where / = y/e. In order to apply theorem [6] we denote it = In r and [/ (u) = 
\Js(e u ). Then taking V = 1 we obtain the generating pair (F, G) for equation 
(JSJ) in the desirable form 

F=U(u), G=^. (12) 

The analytic function $ (from theorem [6]) corresponding to the polar coor- 
dinate system has the form <&(z) =\nz and consequently § z { z ) =1/2. We 
note that $ 2 has a pole in the origin and a zero at infinity. Thus, theorem 
[6] is applicable in any domain Q which does not include these two points. 
Moreover, as for constructing formal powers we need to use the recursive 
integration defined by (TTUT) in what follows we require Q to be any bounded 
simply connected domain not containing the origin. 

From theorem [6] we have that a generating sequence corresponding to the 
generating pair (112"]) can be defined as follows 



U i 
F m = — and G m = — — for even m 

"* £ m z m U 

and 

and G m = — for odd m. 



m z m U m z m 
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As was explained in subsection 13.11 in order to have a complete system of 
formal powers for each n we need to construct Z^ n '(l, z ; z) and Z^ n \i, z ; z). 
For n = we have 

Z^(l,z ;z)=X ( f ) F(z)+^ ) G(z) 

and 

Z^^z ;z) = Xf >) F(z)+^ ) G(z) 
where Aj° , /4°^ are real constants chosen so that 

\f i F{z Q )+$ ) G{z ) = l 
and A,- ' 1 , fif^ are real constants such that 

Xf ) F(z )+^ ) G(z )=z. 
Taking into account that F is real and G is imaginary we obtain that 

A (0) - 1 „(°) - n 

Af = 0, = F(z ). 



Thus, 



and 



F{zq) V e(r ) 



Z (0)( 



zF(zo) . /e(r ) 



z, z ; z) = = i 



F(z) \l e(r) 
where r = \z \. 

For constructing Z^\l, zq; z) and Z^\i, zq; z) we need first the formal 
powers z[°\l, z ; z) and zf\i, z G ; z). According to definition [3] they have 
the form 

z[°\l,z ;z) = X?F 1 (z)+ f i?G 1 (z) 

and 

Z?\i } z,-z) = Xf ) F l (z)+^ ) G 1 (z) 



12 



where A^, ^ are real numbers such that 

X?F 1 (z )+^ ) G 1 (z Q )^l 
and \f\ are real numbers such that 

X^ ) F 1 (z )+^ ) G 1 (z )^i. 

Thus in order to determine \^\ ^ and \ f\ ^ we should solve two systems 
of linear algebraic equations: 

A (D 1 | ..(i) jf^Vo) _i 
1 z e^(r ) + ^ z - L 

and 

A« I + n (1) ' £l/2(ro) = 2 

1 z e^(r ) + ^ z 

which can be rewritten as follows 

A^ +^\e(r )=e 1 / 2 (r )z 

and 

From here we obtain 

\? = e ^(ro)xo, ^ = ^ V VoK A, (1) = -^ 2 (ro)2/o, ^ = ^ 1/2 (r )^ . 
Let us notice that in general for odd m we have 

^£ 0) (i, ^o; *) = -^t^tt + ^ £ (r) , 

7 (o )r . A< m > , vxi"V/*(r) 



» £ l/2( 7 



where 



o- 



A (m) = e l/2( r() ) Rez m = £ l/2( ro ) r m cosm ^ ( 

/iS m) = £- 1/2 (r )Imz m = e- 1/2 (r )rf sin m6 , 
Aj m) = -£ 1/2 (r )Im^ = -e 1/2 {r )r™smm6 , 
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/4 m) = e- 1/2 (r ) Rez™ = £~ 1/2 (r )r™ cosm# 

#o is the argument of the complex number Zq. 
Thus, for odd m: 



Z^(l,2o;,z) = ^— J I cosm8 J , - + isinmfl '' 



z / \ V e(r) V e(r Q ) J ' 



Z<?(i,zb;z) = (-) ( -sinm^oW + i cos m9 J 
In a similar way we obtain the corresponding formulas for even m: 



Z^ } (1,2; ;2;) = J I cos m6» J -^—^ + i sin m6 < 



le{r ) 



e r 



ZW(i,z ;z)= (-j) I -sin m6 J cos m6 J -^y 

In order to apply formula (flOj) for constructing formal powers of higher orders 
we need to calculate the adjoint generating pairs For odd m we 

have 



F* = - G* = e 1/2 (r)z r ' 

m £ i/2( r y m v / 



For even m we obtain 



F* = -tz m e l l 2 (r), G 



' m e i/2( r )- 

Now the whole procedure of construction of formal powers can be easily 
algorithmized. The obtained system of formal powers 



{Z^(l,z ;z), Z^(i,z Q ;z)} 



oo 
n=0 



is complete in the space of all solutions of the main Vekua equation (|S]) with 
/ = £ 1 / 2 (r), i.e., any regular solution W of ([S]) in can be represented in 
the form of a normally convergent series 



oo 



W(Z)=J2 Zin) On, Z ;z) = J2 K Z(n) ^ ^ Z ) + < Z(n) & ^ Z )) 
n=0 n=0 

where a' n = Rea n , a" = Ima n and z is an arbitrary fixed point in Q. 
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